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In most introductory courses on quantum mechanics one is taught that the Hamil- 
tonian operator must be Hermitian in order that the energy levels be real and that 
the theory be unitary (probability conserving). To express the Hermiticity of a 
■ Hamiltonian, one writes H = H\ where the symbol f denotes the usual Dirac Her- 

^ . mitian conjugation; that is, transpose and complex conjugate. In the past few years 

. it has been recognized that the requirement of Hermiticity, which is often stated as 

\ an axiom of quantum mechanics, may be replaced by the less mathematical and more 

physical requirement of space-time reflection symmetry {VT symmetry) without los- 
ing any of the essential physical features of quantum mechanics. Theories defined by 
non-Hermitian T^T-symmetric Hamiltonians exhibit strange and unexpected prop- 
erties at the classical as well as at the quantum level. This paper explains how 
J> \ the requirement of Hermiticity can be evaded and discusses the properties of some 

' non-Hermitian PT-symmetric quantum theories. 

o 

PACS numbers: ll.30.Er, 03.65.-w, 03.65.Bz 

o ' 
o 

The field of PT-symmetric quantum theory is only six years old but already hundreds 
^ ' of papers have been published on various aspects of PT-symmetric quantum mechanics 
and "PT-symmetric quantum field theory. Three international conferences have been held 
Q^. (Prague, 2003; Prague, 2004; Shizuoka, 2004) and three more conferences are planned. Work 

on VT symmetry began with the investigation of quantum-mechanical models and has now 
extended into many areas including quasi-exact solvability, supersymmetry, and quantum 
field theory. Recently, it has been recognized that there is a connection between PT- 
symmetric quantum mechanics and integrable models. The aim of this paper is to introduce 
the subject at an elementary level and to elucidate the properties of theories described by 
PT-symmetric Hamiltonians. This paper will make the field of VT symmetry accessible to 
students who are interested in exploring this exciting, new, and active area of physics. 

The central idea of PT-symmetric quantum theory is to replace the condition that the 
Hamiltonian of a quantum theory be Hermitian with the weaker condition that it possess 
space-time reflection symmetry (PT symmetry). This allows us to construct and study 
many new kinds of Hamiltonians that would previously have been ignored. These new 
Hamiltonians have remarkable mathematical properties and it may well turn out that these 
new Hamiltonians will be useful in describing the physical world. It is crucial, of course, 
that in replacing the condition of Hermiticity by PT symmetry we do not give up any of 
the key physical properties that a quantum theory must have. We will see that if the PT 
symmetry of the Hamiltonian is not broken, then the Hamiltonian will exhibit all of the 
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features of a quantum theory described by a Hermitian Hamiltonian. (The word broken as 
used here is a technical term that will be explained in Sec. ITTl ) 

Let us begin by reviewing some basic ideas of quantum theory. For simplicity, in this 
paper we restrict our attention to one-dimensional quantum-mechanical systems. Also, we 
work in units where Planck's constant h = 1. In elementary courses on quantum mechanics 
one learns that a quantum theory is specified by the Hamiltonian operator that acts on a 
Hilbert space. The Hamiltonian H does three things: 

(i) The Hamiltonian determines the energy eigenstates \En)- These states are the eigen- 
states of the Hamiltonian operator and they solve the time- independent Schrodinger equation 
H\En) = En\En). The energy eigenstates span the Hilbert space of physical state vectors. 
The eigenvalues En are the energy levels of the quantum theory. In principle, one can ob- 
serve or measure these energy levels. The outcome of such a physical measurement is a real 
number, so it is essential that these energy eigenvalues be real. 

(ii) The Hamiltonian H determines the time evolution in the theory. States \t) in the 
Schrodinger picture evolve in time according to the time-dependent Schrodinger equation 
H\t) = —i-^\t), whose formal solution is \t) = e'^^^\0). Operators A(t) in the Heisenberg 
picture evolve according to the time-dependent Schrodinger equation -^A{t) = —i[A{t),H], 
whose formal solution is A(t) = e*^M(0)e~*^*. 

(iii) The Hamiltonian incorporates the symmetries of the theory. A quantum theory 
may have two kinds of symmetries: continuous symmetries, such as Lorentz invariance, and 
discrete symmetries, such as parity invariance and time reversal invariance. A quantum 
theory is symmetric under a transformation represented by an operator A if A commutes 
with the Hamiltonian that describes the quantum theory: [A, H] = 0. Note that if a 
symmetry transformation is represented by a linear operator A and if A commutes with the 
Hamiltonian, then the eigenstates of H are also eigenstates of A. Two important discrete 
symmetry operators are parity (space reflection), which is represented by the symbol V, and 
time reversal, which is represented by the symbol T. The operators V and T are defined 
by their effects on the dynamical variables x (the position operator) and p (the momentum 
operator). The operator V is linear and has the effect of changing the sign of the momentum 
operator p and the position operator f : p — >■ —p and ct — > —x. The operator T is antilinear 
and has the effect p — * —p, x x, and i — > —i. Note that T changes the sign of i because 
(like V) T is required to preserve the fundamental commutation relation [x,p] = z of the 
dynamical variables in quantum mechanics. 

Quantum mechanics is an association between states in a mathematical Hilbert space 
and experimentally measurable probabilities. The norm of a vector in the Hilbert space 
must be positive because this norm is a probability and a probability must be real and 
positive. Furthermore, the inner product between any two different vectors in the Hilbert 
space must be constant in time because probability is conserved. The requirement that the 
probability not change with time is called unitarity. Unitarity is a fundamental property of 
any quantum theory and must not be violated. 

To summarize the discussion so far, the two crucial properties of any quantum theory 
are that the energy levels must be real and that the time evolution must be unitary. There 
is a simple mathematical condition on the Hamiltonian that guarantees the reality of the 
energy eigenvalues and the unitarity of the time evolution; namely, that the Hamiltonian be 
real and symmetric. To explain the term symmetric, as it is used here, let us first consider 
the possibility that the quantum system has only a finite number of states. In this case the 
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Hamiltonian is a finite-dimensional symmetric matrix 

(1) 



H 



(a b c ■ ■ ■ \ 
b d e ■ ■ ■ 
c e f ■■■ 

V; ; ; -.y 



whose entries a, b, c, d, e, f , ■ ■ ■ , are real numbers. For systems having an infinite number 
of states we express H in coordinate space in terms of the dynamical variables x and p. The 
X operator in coordinate space is a real and symmetric diagonal matrix, all of whose entries 
are the real number x. The p operator in coordinate space is imaginary and anti-symmetric 
because P — when it acts to the right but, as we can see using integration by parts, p 
changes sign p = i-^ when it acts to the left. The operator = — ^ is real and symmetric. 
Thus, any Hamiltonian of the form H = p'^ + V{x) when written in coordinate space is 
real and symmetric. (In this paper we use units in which m = \ and we treat x and p as 
dimensionless.) 

However, the condition that H be real and symmetric is not the most general condition 
that guarantees the reality of the energy levels and the unitarity of the time evolution 
because it excludes the possibility that the Hamiltonian matrix might be complex. Indeed, 
there are many physical applications which require that the Hamiltonian be complex. There 
is a more general condition that guarantees spectral reality and unitary time evolution and 
which includes real, symmetric Hamiltonians as a special case. This condition is known 
as Hermiticity. To express the condition that a complex Hamiltonian H is Hermitian we 
write H = W. The symbol f represents Dirac Hermitian conjugation; that is, combined 
transpose and complex conjugation. The condition that H must exhibit Dirac Hermiticity 
is often taught as an axiom of quantum mechanics. The Hamiltonians H — p"^ + p + V{x) 
and H = p"^ + px + xp + V (x) are complex and nonsymmetric but they are Hermitian. 

In this paper we show that while Hermiticity is sufficient to guarantee the two essential 
properties of quantum mechanics, it is not necessary. We describe here an alternative way 
to construct complex Hamiltonians that still guarantees the reality of the eigenvalues and 
the unitarity of time evolution and which also includes real, symmetric Hamiltonians as a 
special case. We will maintain the symmetry of the Hamiltonians in coordinate space, but 
we will allow the matrix elements to become complex in such a way that the condition of 
space-time reflection symmetry {VT symmetry) is preserved. The new kinds of Hamiltonians 
discussed in this paper arc symmetric and have the property that they commute with the VT 
operator: [H, VT] = 0. In analogy with the property of Hermiticity H = H'^ , we will express 
the property that a Hamiltonian is VT symmetric by using the notation H = H^^ . We 
emphasize that our new kinds of complex Hamiltonians are symmetric in coordinate space 
but are not Hermitian in the Dirac sense. To reiterate, acceptable complex Hamiltonians may 
be either Hermitian H = or PT-symmetric H = H^^, but not both. Real symmetric 
Hamiltonians may be both Hermitian and PT-symmetric. 

Using VT symmetry as an alternative condition to Hermiticity, we can construct infinitely 
many new Hamiltonians that would have been rejected in the past because they are not 
Hermitian. An example of such a T'T-symmetric Hamiltonian is 

H^f + ix^. (2) 

We do not regard the condition of Hermiticity as wrong. Rather, the condition of VT 
symmetry offers the possibility of studying new quantum theories that may even describe 
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measurable physical phenomena. Indeed, non-Hermitian PT-symmetric Hamiltonians have 
already been used to describe such phenomena as the ground state of a quantum system of 
hard spheres Reggeon field theory and the Lee- Yang edge singularity Although 
at the time that they were written these papers were criticized for using Hamiltonians that 
were not Hermitian, we now understand that these Hamiltonians have spectral positivity 
and that the associated quantum theories are unitary because these Hamiltonians are VT- 
symmetric. In physics we should keep an open mind regarding the kinds of theories that we 
are willing to consider. Gell-Mann's "totalitarian principle" states that among the possible 
physical theories "Everything which is not forbidden is compulsory." 

This paper is organized as follows: I discuss in a personal way my discovery of VT- 
symmetric quantum mechanics and give a brief history of the early days of this subject in 
Sec. n] Section IIIII explains how to calculate the energy levels of a "PT-symmetric Hamil- 
tonian. Section HVl describes the classical mechanics of PT-symmetric Hamiltonians. Next, 
in Sec. |V]we show that a Hamiltonian having an unbroken PT symmetry defines a unitary 
quantum theory. The demonstration of unitarity is based on showing that PT-symmetric 
Hamiltonians that have an unbroken PT symmetry also possess a new parity-like symmetry; 
this symmetry is represented by a new operator that we call C. We give in Sec. IVIl a simple 
2x2 matrix illustration of the procedures used in Sec. El In Sec. IVIII we discuss the nature 
of observables in PT-symmetric quantum-mechanical theories. We show how to calculate 
the C operator in Sec. IVIIII In Sec. IIXI we explain why one may regard PT-symmetric 
quantum mechanics as a complex version of ordinary quantum mechanics. Finally, in Sec. |3 
we discuss some possible physical applications of PT-symmetric quantum mechanics. 



II. A PERSONAL HISTORY OF PT SYMMETRY 

My first encounter with a non-Hermitian complex Hamiltonian dates back to the summer 
of 1993. In the course of a private conversation with D. Bessis at GEN Saclay, I learned 
that he and J. Zinn- Justin had noticed that the eigenvalues of the Hamiltonian operator in 
(j2)) seemed to be real and they wondered if the spectrum (the set of energy eigenvalues of 
the Hamiltonian) might be entirely real. (Their interest in the Hamiltonian (j2)) was inspired 
by early work on the Lee- Yang edge singularity Q.) At the time I did not plan to pursue 
this conjecture further because it seemed absurd that a complex non-Hermitian Hamiltonian 
might have real energy levels. 

I did not know it at the time, but Bessis and Zinn- Justin were not the first to notice that a 
complex cubic quantum-mechanical Hamiltonian might have real eigenvalues. For example, 
early studies of Reggeon field theory in the late 1970's led a number of investigators to 
observe that model cubic quantum-mechanical Hamiltonians like that in (j21) might have 
real eigenvalues Q. Also, E. Caliceti et al. observed in 1980 that on the basis of Borel 
summability arguments the spectrum of a Hamiltonian related to (0) is real 0]. In each 
of these cases the possibility that a complex non-Hermitian Hamiltonian might have real 
energy levels was viewed as an isolated curiosity. It was believed that such a Hamiltonian 
could not describe a valid theory of quantum mechanics because the non-Hermiticity of the 
Hamiltonian would result in nonunitary time evolution 

I did not forget about the remarkable Hamiltonian in (j2)) and in 1997 I decided to investi- 
gate it. I suspected that if the spectrum of this Hamiltonian was real, it was probably due to 
the presence of a symmetry and I realized that (j21) does possess PT symmetry because any 
real function of ix is PT-symmetric. I decided that a simple and natural way to determine 
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the spectrum of ^ would be to use the delta expansion, a perturbative technique that I 
had developed several years earlier for solving nonlinear problems I asked my former 
graduate student S. Boettcher to join me in this investigation. 

The delta expansion is an extremely simple technique for solving nonlinear problems 
perturbatively (approximately). The idea of the delta expansion is to introduce a small 
perturbation parameter 6 into a nonlinear problem in such a way that 5 is a measure of 
the nonlinearity of the problem. To illustrate how the delta expansion is used, consider the 
Thomas-Fermi equation, a difficult nonlinear boundary-value problem that describes the 
approximate electric charge distribution in an atom: 

y"{x)=y'/'x-'/', y(0) = l, y(oo) = 0. (3) 

There is no exact closed-form solution to this problem. However, a nice way to solve this 
problem perturbatively is to introduce the parameter S in the exponent: 

y"{x)=y{y/xY, y{0) = l, y(oo) = 0. (4) 

When 5 = |, (0} reduces to (jH))- However, in (j3)) we treat the parameter 6 as small (5 ^ 1). 
When (5 = 0, the problem becomes linear and therefore it can be solved exactly; the solution 
is yo{x) = e~^. This is the leading term in the perturbation expansion for y{x), which has 
the form y{x) = Xl^^o ^"(■^)'^"- ^^^^ ^'^ calculate the coefficients of the higher powers 
of S. At the end of the calculation one sets 6 = 1/2, and from just the first few terms in 
the perturbation series one obtains a good numerical approximation to the solution to the 
Thomas- Fermi equation. 

I was eager to find out what would happen if we applied delta-expansion methods to 0. 
We replaced the Hamiltonian by the one-parameter family of Hamiltonians 

H = p^ + x\txY, (5) 

where 6 is regarded as a small real parameter. There are two advantages in inserting 6 in 
this fashion: First, the new Hamiltonian remains PT-symmetric for all real 6. Thus, the 
insertion of 6 maintains the VT symmetry of the original problem. Second, when 6 = 0, 
the Hamiltonian (0) reduces to that of the harmonic oscillator, which can be solved exactly 
because the underlying classical equations of motion are linear. Each of the energy levels of 
(0) have a delta expansion of the general form 

E = EZo^nS-. (6) 

The series coefficients are easy to calculate and they are real. Thus, assuming that this delta 
expansion converges, the eigenvalues of if in (jS)) must be real. At the end of the calculation 
we set 5 = 1 in © to recover the eigenvalues of the original Hamiltonian in 0. 

The problem with the delta expansion is that it is difficult to prove rigorously that the 
expansion converges. We were able to conclude only that for every n there is always a 
neighborhood about 5 = in which the delta expansion for the first n eigenvalues of (0) 
converges and thus these eigenvalues are real when 6 is real. Our discovery that the first n 
eigenvalues of the complex Hamiltonian (0) were real for a small range of real S near 6 = 
was astonishing to us, but we were disappointed that the delta expansion was not powerful 
enough to determine whether all of the eigenvalues of the original complex Hamiltonian (j21) 
are real. However, our delta-expansion analysis inspired us to perform detailed perturbative 
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and numerical studies of the spectrum of if in (P). To our amazement we found that all of 
the eigenvalues of H remain real for all 6 > [6J . We coined the term PT-symmetric to 
describe these new non-Hermitian complex Hamiltonians having real energy levels 0]. 
To present the results of our numerical studies, we rewrite the Hamiltonian (0) as 

H = p'- {^xf, (7) 

where is a continuous real parameter . The eigenvalues of this Hamiltonian are entirely 
real for all > 2, while for N < 2 the spectrum is partly real and partly complex. Clearly, 
the Hamiltonian H = -\- ix^ in ^ is just one member of a huge and remarkable class of 
complex Hamiltonians whose energy levels are real and positive. The spectrum of H exhibits 
three distinct behaviors as a function of (see Fig. Q): 

(i) When N > 2 the spectrum is infinite, discrete, and entirely real and positive. This 
region includes the case N = 4 for which H = p'^ — x^. Amazingly, the spectrum of this 
wrong-sign potential is positive and discrete. (Also, {x) 7^ in the ground state because 
H breaks parity symmetry!) At the lower bound A^ = 2 of this region lies the harmonic 
oscillator, (ii) A transition occurs at A^ = 2. When 1 < A^ < 2 there are only a finite 
number of positive real eigenvalues and an infinite number of complex conjugate pairs of 
eigenvalues. We say that in this region the VT symmetry is broken and that A^ > 2 is 
a region of unbroken VT symmetry. (We explain the notion of broken and unbroken VT 
symmetry in greater detail below.) As A^ decreases from 2 to 1, adjacent energy levels merge 
into complex conjugate pairs beginning at the high end of the spectrum. Ultimately, the 
only remaining real eigenvalue is the ground-state energy, which diverges as A^ — l"*" jof. 
(iii) When A^ < 1 there are no real eigenvalues. 

It is apparent that the reality of the spectrum of ((Tj) when A^ > 2 is connected with its 
VT symmetry. The association between VT symmetry and the reality of spectra can be 
understood as follows: We say that the VT symmetry of a Hamiltonian H is unbroken if all 



19 

17 

15 

13 

M 11 
u 

7 
5 
3 
1 





/ / / ^ 






- \ 




V 





FIG. 1: Real energy levels of the Hamiltonian H = fp' — {ix)^ as a function of the parameter A^. 
When A^ > 2, the entire spectrum is real and positive. The lower bound of this region, N = 2, 
corresponds to the harmonic oscillator, whose energy levels are En = 2n + 1. When 1 < A^ < 2, 
there are a finite number of positive real eigenvalues and infinitely many complex conjugate pairs 
of eigenvalues. As decreases from 2 to 1, the number of real eigenvalues decreases and when 
A^ < 1.42207, the only real eigenvalue is the ground-state energy. As A^ — > 1"^, the ground-state 
energy diverges. For < 1 there are no real eigenvalues. 



7 



of the eigenfunctions of H are simultaneously eigenfunctions of VT [K 

Here is a proof that if the VT symmetry of a Hamiltonian H is unbroken, then the 
spectrum of H is real: Assume that a Hamiltonian H possesses VT symmetry (i.e., that 
H commutes with the VT operator) and that if is an eigenstate of H with eigenvalue E, 
then it is simultaneously an eigenstate of VT with eigenvalue A: 

H(f) = E(f) and = A0. (8) 

We begin by showing that the eigenvalue A is a pure phase. Multiply VT(j) = Xcf) on the 
left by VT and use the fact that V and T commute and that = = 1 to conclude 
that (f) = X*X(f) and thus A = e*" for some real a. Next, introduce a convention that we use 
throughout this paper. Without loss of generality we replace the eigenfunction by e~*"/^0 
so that its eigenvalue under the operator VT is unity: 

PT0 = 0. (9) 

Next, we multiply the eigenvalue equation H(f) = Ecj) on the left by VT and use \VT , if] = 
to obtain Ecf) = E*(j). Hence, E = E* and the eigenvalue E is real. 

The crucial assumption in this argument is that is simultaneously an eigenstate of H 
and of VT. In quantum mechanics if a linear operator X commutes with the Hamiltonian 
H, then the eigenstates of H are also eigenstates of X. However, we emphasize that the 
operator VT is not linear (it is antilinear) and thus we must make the extra assumption 
that the VT symmetry of H is unbroken; that is, is simultaneously an eigenstate of H and 
VT. This extra assumption is nontrivial because it is hard to determine a priori whether 
the VT symmetry of a particular Hamiltonian H is broken or unbroken. For if in ((7j) the 
VT symmetry is unbroken when N > 2 and it is broken when N < 2. The conventional 
Hermitian Hamiltonian for the quantum-mechanical harmonic oscillator lies at the boundary 
of the unbroken and the broken regimes. 

I am delighted at the research activity that my work has inspired. In 2001 Dorey et 
al. proved rigorously that the spectrum of ii in (|7j) is real and positive in the region > 
2. Dorey et al. used techniques such as the Bethe ansatz and the Baxter-TQ relation, which 
are used in the study of integrable models and conformal quantum field theory. In doing 
so they have helped to establish a remarkable connection between the ordinary differential 
equation (the Schrodinger equation) that describes PT-symmetric quantum mechanics and 
the study of integrable models. This connection, which has become known as the ODE/IM 
correspondence, is rich and profound and will lead to a much deeper understanding of both 
types of theories. Many other PT-symmetric Hamiltonians for which space-time reflection 
symmetry is not broken have been investigated, and the spectra of these Hamiltonians have 
also been shown to be real and positive 0, 0, 0, 3|. Evidently, the phenomenon of VT 



symmetry is quite widespread and arises in many contexts. 



III. ENERGY LEVELS OF A PT- SYMMETRIC HAMILTONIAN 

The purpose of this section is to explain how to calculate the eigenvalues of the complex 
Hamiltonian operator in (jHI). To calculate the energy levels of a PT-symmetric Hamilto- 
nian we adopt the techniques that are used for calculating the energy levels of conventional 
Hermitian Hamiltonians. These techniques involve converting the formal eigenvalue prob- 
lem ii0 = Ecj) to a Schrodinger differential equation whose solutions satisfy appropriate 
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FIG. 2: Wedges in the complex-x plane containing the contour on which the eigenvalue problem 
for the differential equation H10|) for N = 4.2 is posed. In these wedges (j){x) vanishes exponentially 
as oo. The wedges are bounded by lines along which the solution to the differential equation 

is oscillatory. 



boundary conditions. This Schrodinger equation is then solved numerically or by using 
approximate methods such as WKB. 

The Schrodinger eigenvalue problem for the PT-symmetric Hamiltonian ((Tj) is 



- 0AX - [IX 



(10) 



where En is the nth eigenvalue. For a Hermitian Hamiltonian the boundary conditions 



as 



x\ 



oo on 



that give quantized energy levels En are that the eigenfunctions 
the real axis. This condition suffices for (|l(Jj) when 1 < < 4, but for > 4 we must 
continue the eigenvalue problem for (fTIH) into the complex-x plane. Thus, we replace the 
real-x axis by a contour in the complex plane along which the differential equation holds. 
The boundary conditions that lead to quantization are imposed at the endpoints of this 
contour. (The rest of this brief section is somewhat technical and may be skipped by those 
without a background in complex differential equations. See Ref. for more information 
on how to solve such problems.) 

The endpoints of this contour lie in regions in the complex-x plane in which the eigen- 
functions (pnix) —>■ exponentially as |a;| — *• oo. These regions are known as Stokes wedges 
(see Fig. The Stokes wedges are bounded by the lines along which the solution to the 
differential equation is oscillatory There are many wedges in which we can require 

that 0(x) — as |x| oo. Thus, there are many eigenvalue problems associated with a 
given differential equation jl^. For a given value of we must first identify which one 
of these eigenvalue problems is associated with (fTUI) . To do so we start with the harmonic 
oscillator problem at = 2 and smoothly vary the parameter A^ until it reaches the given 
value. At A^ = 2 the eigenfunctions vanish in wedges of angular opening centered about 
the negative-real and positive-real x axes. For any A^ > 1 the centers of the left and right 
wedges lie at the angles 
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left 



-TT 



N-2 

2N + 4: 



TC 



and 9ri 



right 



jV-2 

' 2N+4: 



TT. 



:ii) 



The opening angle of these wedges is A 



N+2 



IT. 



The differential equation (jlUj) may be 



integrated on any path in the complex-x plane so long as the path approaches complex 
infinity inside the left wedge and inside the right wedge. These wedges contain the real- a; 
axis when 1 < A^ < 4. 
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As increases from 2, the left and right wedges rotate downward into the complex-x 
plane and become thinner. We can see on Fig.^that the eigenvalues grow with N as N ^ oo. 
At = oo the differential equation contour runs up and down the negative imaginary axis, 
and this leads to an interesting limiting eigenvalue problem. Because all of the eigenvalues 
of ((Zj) diverge like A^^ as A^ ^ cxd, for large A^ we replace H by the rescaled Hamiltonian 
H/N"^. In the limit N —>■ oo this new Hamiltonian becomes exactly solvable in terms of 
Bessel functions. The eigenvalue problem at A^ = cxd is the PT-symmetric equivalent of the 
square well in ordinary Hermitian quantum mechanics [l6| . 

As A^ decreases below 2, the wedges become wider and rotate into the upper-half x plane. 
At A^ = 1 the angular opening of the wedges is |7r and the wedges are centered at |7r and ^vr. 
Thus, the wedges become contiguous at the positive-imaginary x axis, and the differential 
equation contour can be pushed off to infinity. Hence, there is no eigenvalue problem when 
A^ = 1 and, as we would expect, the ground-state energy diverges as A^ — > 1"^ (see Fig. [T)). 

Having defined the eigenvalue problem for the Hamiltonian in ((7j), we can solve the dif- 
ferential equation by using numerical methods. We can also can use approximate analytical 
methods such as WKB (l rj. WKB gives a good approximation to the eigenvalues in Fig. ^ 
when N > 2. The novelty of this WKB calculation is that it must be performed in the 
complex plane. The turning points x± are those roots of E + {ix)^ = that analytically 
continue off the real A^ moves away from N = 2: 

^_^^l/iVgM3/2-l/iV)^ ^^^^l/iVg-..(l/2-l/iV)_ (^2) 

These points lie in the lower (upper) x plane in Fig. |21when N > 2 {N < 2). 

The WKB quantization condition is {n + l/2)7r = J^'^ dx \/ E + (ix)^ . It is crucial that 
the integration path be such that this integral is real. When N > 2 this path lies entirely 
in the lower-half x plane, and when N = 2 the path lies on the real axis. When N < 2 the 
path is in the upper-half x plane; it crosses the cut on the positive-imaginary axis and thus 
is not a continuous path joining the turning points. Hence, WKB fails when N < 2. 

When N >2, the WKB calculation gives 



En ~ 



r(3/2 + l/Ar)0F(n + l/2) ]^^/(^+^) 

sin(vr/Ar)r(l + l/AT) J ^ ^^^^ ^'^^ 



This result is quite accurate. The fourth exact eigenvalue (obtained using Runge-Kutta) for 
the case A^ = 3 is 11.3143 while WKB gives 11.3042, and the fourth exact eigenvalue for the 
case N = Ais 18.4590 while WKB gives 18.4321. 



IV. T'T-SYMMETRIC CLASSICAL MECHANICS 

In the study of classical mechanics the objective is to describe the motion of a particle 
satisfying Newton's law F = ma. The trajectory x{t) of the particle is a real function of 
time t. The classical equation of motion for the complex PT-symmetric Hamiltonian (|7j) 
describes a particle of energy E subject to complex forces. Thus, we have the surprising 
result that classical "PT-symmetric Hamiltonians describe motion that is not limited to the 
real-x axis. The classical path x{t) may lie in the complex-x plane. The purpose of this 
section is to describe this remarkable possibility fis']. 

An intriguing aspect of Fig. ^ is the transition at A^ = 2. As A^ goes below 2, the 
eigenvalues begin to merge into complex conjugate pairs. The onset of eigenvalue merging 
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FIG. 3: Classical paths in the complex-x plane for the N = 2 oscillator. The paths form a set of 
nested ellipses. These closed periodic orbits occur when the VT symmetry is unbroken. 
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FIG. 4: Classical paths in the complex-x plane for the = 3 oscillator. In addition to the periodic 
orbits, one path runs off to ioo from the turning point on the imaginary axis. 



can be thought of as a phase transition. We show in this section that the underlying cause 
of this quantum transition can be understood by studying the theory at classical level. 
The trajectory x{t) of a classical particle governed by the PT-symmetric Hamiltonian 



dZj) obeys ±dx[E + 



IX 



,Afl-l/2 



2dt. While E and dt are real, x{t) lies in the complex plane 



in Fig. 13 When N = 2 (the harmonic oscillator), there is one classical path that terminates 
at the classical turning points x± in (fT2|) . Other paths are nested ellipses with foci at the 
turning points (see Fig. |2I). All these paths have the same period. 

When = 3, there is again a classical path that joins the left and right turning points 
and an infinite class of paths enclosing the turning points (see Fig. H}. As these paths 
increase in size, they approach a cardioid shape (see Fig.EJ. The indentation in the limiting 
cardioid occurs because paths may not cross, and thus all periodic paths must avoid the 
path in Fig. |3] that runs up the imaginary axis. When is noninteger, we obtain classical 
paths that move off onto different sheets of the Riemann surface (see Fig. IH)). 

In general, whenever N >2, the trajectory joining x± is a smile-shaped arc in the lower 
complex plane. The motion is periodic. Thus, the equation describes a complex pendulum 
whose (real) period T is given by 



T 



2-JV 

2E COS 



{N-2)n 
2N 



T{l/2 + l/N) ■ 



(14) 



Below the transition at = 2 a path starting at one turning point, say x_, moves toward 
but misses the turning point x+. This path spirals outward, crossing from sheet to sheet 
on the Riemann surface, and eventually veers off to infinity. Hence, the period abruptly 
becomes infinite. The total angular rotation of the spiral is finite for all A^ < 2, but it 
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FIG. 5: Classical paths in the complex-x plane for the = 3 oscillator. As the paths get larger, 
they approach a shape resembling a cardioid. We have plotted the resettled paths. 

becomes infinite as — 2~ (see Fig. Ej). 

V. PT- SYMMETRIC QUANTUM MECHANICS 

The discovery that the eigenvalues of many PT-symmetric Hamiltonians are real and 
positive raises an urgent question: Does a non-Hermitian Hamiltonian such as if in ((Zj) 
define a physical theory of quantum mechanics or is the positivity of the spectrum merely 
an intriguing mathematical property of special classes of complex eigenvalue problems? A 
physical quantum theory (i) must possess a Hilbert space of state vectors and this Hilbert 
space must have an inner product with a positive norm; (ii) the time evolution of the theory 
must be unitary; that is, the norm must be preserved in time. 

A definitive answer to this question has been found |[19|,(20|]. For a complex non-Hermitian 
Hamiltonian having an unbroken VT symmetry, a linear operator C that commutes with 
both H and VT can be constructed. We denote the operator representing this symmetry 
by C because the properties of C are similar to those of the charge conjugation operator in 
ordinary particle physics. The inner product with respect to CVT conjugation is 

(V'lx)'^^ = /^a;^^^^(x)x(x), (15) 

where ip'^'^'^{x) = J dyC{x,y)ip*{—y). This inner product satisfies the requirements for the 
quantum theory defined by H to have a Hilbert space with a positive norm and to be a 
unitary theory of quantum mechanics. 

To explain the construction of the C operator we begin by summarizing the mathematical 
properties of the solution to the eigenvalue problem (fTUj) associated with the Hamiltonian H 
in d?!). Recall from Fig. |21 that this differential equation is imposed on an infinite contour in 
the complex-x plane and that for large |x| this contour lies in wedges placed symmetrically 
with respect to the imaginary-x axis as in Fig. |21 When N > 2, H has an unbroken 
VT symmetry. Thus, the eigenfunctions (f)n{x) are simultaneously eigenstates of the VT 
operator: VT(j)n{x) = A„(^„(x). As we argued in Sec. HIl A„ can be absorbed into 0„(x) so 
that VT (j)^{x) = (pli-x) = 0„(x) [see ©]. 

The eigenstates of a conventional Hermitian Hamiltonian are complete. There is strong 
evidence that the eigenfunctions (pnix) for the PT-symmetric Hamiltonian (|7j) are also com- 
plete. The coordinate-space statement of completeness is 



E~=o(-l)"'^"(^)<^"(l/) =6{x-y) (x, y real). 



(16) 
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FIG. 6: Classical paths for the case N = 2.5. These paths do not intersect. The graph shows the 
projection of the parts of the path that lie on three different sheets of the Riemann surface. As the 
size of the paths increases, a limiting cardioid appears on the principal sheet. On the remaining 
sheets of the surface the path exhibits a knot-like topological structure. 



This nontrivial result has been verified numerically to extremely high accuracy (twenty deci- 
mal places) 2^1. The unusual factor of (—1)" in this sum does not appear in conventional 



quantum mechanics. This factor is explained in the following discussion. 

We must now try to find the inner product associated with our PT-symmetric Hamil- 
tonian and it is here that we can see the difficulty connected with its non-Hermiticity. In 
conventional Hermitian quantum mechanics the Hilbert space inner product is specified 
even before we begin to look for the eigenstates of H. For our non-Hermitian Hamiltonian 
we must try to guess the inner product. A reasonable guess for the inner product of two 
functions f{x) and g{x) might be 

{f,g)^Jdx[VTf{xMx), (17) 

where VTf{x) = [f{—x)]* and the path of integration in the complex-x plane follows the 
contour described in Sec. IIIII The apparent advantage of this choice for the inner product 
is that the associated VT norm (/, /) is independent of the overall phase of f{x) and is 
conserved in time. With respect to this inner product the eigenfunctions 4>m{x) and 0„(a;) 
of in ([7j) are orthogonal for n ^ m. However, when m = n we see that the VT norms of 
the eigenfunctions are not positive: 

(0m,0n) = (-l)"5™„. (18) 

This result is apparently true for all values of in ([7j) and it has been verified numerically 
to extremely high precision. Because the norms of the eigenfunctions alternate in sign, the 
metric associated with the VT inner product (■, ■) is indefinite. This sign alternation is a 
generic feature of the VT inner product. Extensive numerical calculations verify that (jl8p 
holds for all > 2. 

Despite the existence of a nonpositive inner product, we can still do some of the analysis 
that one would normally perform for a conventional Schrodinger equation H(j)n = En4>n- For 
example, we can use the inner product formula ()18|) to verify that ()16p is the representation 
of the unity operator by showing that / dy6{x — y)S{y — z) = 6{x — z). We can use 
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FIG. 7: Classical paths in the complex-x plane for N = 1.8, N = 1.85 and N = 1.9. These 
nonperiodic paths spiral outward to infinity. As ^ 2 from below, the number of turns in the 
spiral increases. The lack of periodic orbits corresponds to a broken VT symmetry. 

completeness to reconstruct the parity operator V in terms of the eigenstates. The parity 
operator in position space is 

n^, y) = En=oi-'^rM^)M-y) = ^(^ + y)- (i9) 

By virtue of (fTHj) the square of the parity operator is unity: V"^ = 1. We can also reconstruct 
H in coordinate space: H{x, y) = 1)'^-E„0„(a:;)0„(?/). Using (fTBj) - (fTHj) we can see that 

this Hamiltonian satisfies H(j)n{x) = En(f)n{x)- 

We now address the question of whether a PT-symmetric Hamiltonian defines a physically 
viable quantum mechanics. The difficulty with formulating a PT-symmetric quantum theory 
is that the vector space of quantum states is spanned by the energy eigenstates, of which 
half have norm +1 and half have norm —1. In quantum theory the norms of states carry a 
probabilistic interpretation, so the indefinite metric (jl8p is unacceptable. 

The situation here in which half of the energy eigenstates have positive norm and half 
have negative norm is analogous to the problem that Dirac encountered in formulating the 
spinor wave equation in relativistic quantum theory [2^1 . Following Dirac, we attack the 
problem of an indefinite norm by finding an interpretation of the negative-norm states. We 
claim that in any theory having an unbroken VT symmetry there exists a symmetry of the 
Hamiltonian connected with the fact that there are equal numbers of positive- and negative- 
norm states. To describe this symmetry we construct the aforementioned linear operator C 
in position space as a sum over the eigenstates of the Hamiltonian jl9l |: 

C{x,y) = j:n=oU^)Uy)- (20) 

The properties of this new operator C resemble those of the charge conjugation operator 
in quantum field theory. For example, we can use (fTH|) - (fTH|) to verify that the square of C 
is unity (C^ = 1): / dyC{x, y)C{y, z) = 6{x — z). Thus, the eigenvalues of C are ±1. Also, C 
commutes with the Hamiltonian H. Therefore, since C is linear, the eigenstates of H have 
definite values of C. Specifically, if the energy eigenstates satisfy (fTHj) . then we have 

C(f)n{x) = JdyC{x,y)(j)n{y) = TZ=o(t>m{x) j dy (l):m{y)(pn{y) = (-l)"0„(x). 
Thus, C represents the measurement of the sign of the VT norm of an eigenstate in (|TH)). 
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The operators V and C are distinct square roots of the unity operator 5{x — y). That is, 
'P'^ = = 1^ but V ^ C. Indeed, V is real, while C is complex. The parity operator in 
coordinate space is explicitly real V{x,y) = 6{x + y), while the operator C{x,y) is complex 
because it is a sum of products of complex functions, as we see in (QUI) . The two operators 
V and C do not commute. However, C does commute with VT. 

Finally, having obtained the operator C we define the new inner product structure given 
in fll5j) . This inner product has a positive definite norm. Like the VT inner product ()17|1 
this new inner product is phase independent. Also, it is conserved in time because the 
time evolution operator (just as in ordinary quantum mechanics) is e*^*. The fact that 
H commutes with VT and with CVT implies that both inner products, p7|) and (fT3jl . 
remain time independent as the states evolve. However, unlike fll7|) . the inner product (fT^ 
is positive definite because C contributes —1 when it acts on states with negative VT norm. 
In terms of the CVT conjugate, the completeness condition ([THjl reads 

En=oU^)[C'P'^My)] = ^(^ - y)- (21) 

To review, in the mathematical formulation of a conventional quantum theory the Hilbert 
space of physical states is specified first. The inner product in this vector space is defined 
with respect to ordinary Dirac Hermitian conjugation (complex conjugate and transpose). 
The Hamiltonian is then chosen and the eigenvectors and eigenvalues of the Hamiltonian are 
determined. In contrast, the inner product for a quantum theory defined by a non-Hermitian 
PT-symmetric Hamiltonian depends on the Hamiltonian itself and is thus determined dy- 
namically. One can view this new kind of quantum theory as a "bootstrap" theory because 
one must solve for the eigenstates of H before knowing what the Hilbert space and the 
associated inner product of the theory are. The Hilbert space and the CVT inner product 
(fT3|) are then determined by these eigenstates via (pUj). 

The operator C does not exist as a distinct entity in ordinary Hermitian quantum me- 
chanics. Indeed, if we allow the parameter iV in ((7j) to tend to 2, the operator C in this limit 
becomes identical to V. Thus, in this limit the CVT operator becomes T, which is just 
complex conjugation. As a consequence, the inner product (fT3|) defined with respect to the 
CVT conjugation reduces to the complex conjugate inner product of conventional quantum 
mechanics when N 2. Similarly, in this limit ()21|) reduces to the usual statement of 
completeness ^n^n{x)4>*n{y) = 6{x -y). 

The CVT inner-product (fT^ is independent of the choice of integration contour C as 
long as C lies inside the asymptotic wedges associated with the boundary conditions for 
the eigenvalue problem (jlUj) . In ordinary quantum mechanics, where the positive-definite 
inner product has the form J dx f*{x)g{x), the integral must be taken along the real axis 
and the path of the inte grat ion cannot be deformed into the complex plane because the 
integrand is not analytic j25|. The VT inner product ()17p shares with ()15p the advantage 
of analyticity and path independence, but it suffers from nonpositivity. It is surprising that 
we can construct a positive-definite metric by using CVT conjugation without disturbing 
the path independence of the inner-product integral. 

We can now explain why PT-symmetric theories are unitary. Time evolution is expressed 
by the operator e"*"^*, whether the theory is determined by a PT-symmetric Hamiltonian 
or just an ordinary Hermitian Hamiltonian. To establish unitarity we must show that as 
a state vector evolves, its norm does not change in time. If ipo{x) is a prescribed initial 
wave function belonging to the Hilbert space spanned by the energy eigenstates, then it 
evolves into the state iptix) at time t according to iptix) = e~^^^ipQ{x). With respect to the 
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CVT inner product defined in ()15|) . the norm of tlie vector iptix) does not cliange in time, 
(iptl'ipt) = (V'olV'o)) because tlie Hamiltonian H commutes witli the CVT operator. 

VI. ILLUSTRATIVE EXAMPLE: A 2 x 2 MATRIX HAMILTONIAN 

The 2x2 matrix Hamiltonian 

^=(''f.e--0 (22) 

where the three parameters r, s, and 6 are real, illustrates the above results on PT-symmetric 
quantum mechanics. This Hamiltonian is not Hermitian, but it is VT symmetric, where the 
parity operator is V = (° o) ^ performs complex conjugation (26l |. 

There are two parametric regions for this Hamiltonian. When < sin^ 9, the energy 
eigenvalues form a complex conjugate pair. This is the region of broken VT symmetry. On 
the other hand, when > sin^ 9, then the eigenvalues e± = r cos 6' ± (s^ — sin^ 9Y^'^ are 
real. This is the region of unbroken VT symmetry. In the unbroken region the simultaneous 
eigenstates of the operators H and VT are 

= ( ) and |g-)= ' (23) 

V 2 cos a \ ^ / V 2 cos a \ ^ J 

where we set sin a = (r/s) sin 6*. The VT inner product gives (e±,e±) = ±1 and {e±^e^) = 
0, where {u,v) = (VTu) - v. Therefore, with respect to the VT inner product, the resulting 
vector space spanned by the energy eigenstates has a metric of signature (+,—). The 
condition > sin^ 9 ensures that the VT symmetry is not broken. If this condition 
is violated, the states are no longer eigenstates of VT because a becomes imaginary. 
When VT symmetry is broken, the VT norm of the energy eigenstate vanishes. 
Next, we construct the operator C using (j^ : 

^ _ 1 / i sin a 1 
cos a V -^sina 

Note that C is distinct from H and V and it has the key property that C\e±) = ±\6±). The 
operator C commutes with H and satisfies = 1. The eigenvalues of C are precisely the 
signs of the VT norms of the corresponding eigenstates. Using the operator C we construct 
the new inner product structure {u\v) = (CVTu) ■ v. This inner product is positive definite 
because {e±\e±) = 1. Thus, the two-dimensional Hilbert space spanned by \e±), with inner 
product (•!•), has signature (+,+). 

Finally, we show that the CVT norm of any vector is positive. For the arbitrary vector 
ip = (^), where a and b are any complex numbers, we see that 

\o* / \0' J COS a \o* — ta* sm a / 

Thus, (iplip) = {CVTip) ■ ip = ■^^[ci*ci + b*b + i{b*b — a* a) sin a]. Now \et a = x + iy and 
b = u + iv, where x, y, u, and v are real. Then 




= (x^ + f ^ + 2xv sin a + y'^ + — 2yu sin a) , 

cos a ^ ^ 



(25) 
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which is exphcitly positive and vanishes only ifx = y = u = v = 0. 

Since {u\ denotes the CPT-conjugate of \u), the completeness condition reads 

k+)(^+l + k-)(^-l = (J (26) 

Furthermore, using the CVT conjugate (£±1, we get C as C = — \e-){e-\. 

If we set 6' = in this two-state system, the Hamiltonian becomes Hermitian. 
However, C then reduces to the parity operator V. As a consequence, CVT invariance 
reduces to the standard condition of Hermiticity for a symmetric matrix; namely, that 
H = H*. This is why the hidden symmetry C was not noticed previously. The operator C 
emerges only when we extend a real symmetric Hamiltonian into the complex domain. 

VII. OBSERVABLES IN PT-SYMMETRIC QUANTUM MECHANICS 

How do we represent an observable in PT-symmetric quantum mechanics? Recall that in 
ordinary quantum mechanics the condition for a linear operator A to be an observable is that 
A = . This condition guarantees that the expectation value of A in a state is real. Because 
operators in the Heisenberg picture evolve in time according to A{t) = e'^*yl(0)e~*^*, this 
Hermiticity condition is maintained in time. In PT-symmetric quantum mechanics the 
equivalent condition is that at time t = the operator A must obey the condition A"^ = 
CVTACVT, where A"^ is the transpose of A. If this condition holds at t = 0, then it will 
continue to hold for all time because we have assumed that H is symmetric {H = H'^). This 
condition also guarantees that the expectation value of A in any state is real. 

The operator C itself satisfies this requirement, so it is an observable. Also the Hamil- 
tonian is an observable. However, the x and p operators are not observables. Indeed, the 
expectation value of x in the ground state is a negative imaginary number. Thus, there is 
no position operator in PT-symmetric quantum mechanics. In this sense PT-symmetric 
quantum mechanics is similar to fermionic quantum field theories. In such theories the 
fermion field corresponds to the x operator. The fermion field is complex and does not have 
a classical limit. One cannot measure the position of an electron; one can only measure the 
position of the charge or of the energy of the electron! 

One can see why the expectation of the x operator is a negative imaginary number by 
examining Fig. Note that the classical trajectories have left-right {VT) symmetry, but 
not up-down symmetry. Also, the classical paths favor the lower-half complex-x plane. 
Thus, the average classical position is a negative imaginary number. Just as the classical 
particle moves about in the complex plane, the quantum probability current flows about 
in the complex plane. It may be that the correct interpretation is to view PT-symmetric 
quantum mechanics as describing the interaction of extended, rather than pointlike objects. 

VIII. CALCULATION OF THE C OPERATOR 

The distinguishing feature of PT-symmetric quantum mechanics is the C operator. The 
discovery of C in Ref. [l^ raises the question of how to evaluate the formal sum in fOUl) that 
represents C. In ordinary Hermitian quantum mechanics there is no such operator. Only a 
non- Hermitian PT-symmetric Hamiltonian possesses a C operator distinct from the parity 
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operator V. Indeed, if we were to evaluate ()2U|) for a PT-symmetric Hamiltonian that is 
also Hermitian, the result would be V, which in coordinate space is 5{x + y) [see (UHl)]- 

Calculating C by direct evaluation of the sum in ()20j) is not easy in quantum mechanics 
because it is necessary to calculate all the eigenfunctions of H. Such a procedure 

cannot be used in quantum field theory because there is no simple analog of the Schrodinger 
eigenvalue differential equation and its associated coordinate-space eigenfunctions. 

Fortunately, there is an easy way to calculate the C operator, and the procedure circum- 
vents the difficult problem of evaluating the sum in ()2U|). As a result the technique readily 
generalizes from quantum mechanics to quantum field theory. In this section we use this 
technique to calculate C for the "PT-symmetric Hamiltonian [27[ 

H = \p^ + + lex^. (27) 

We will show how to calculate C perturbatively to high order in powers of e for this cubic 
Hamiltonian. Calculating C for other kinds of interactions is a bit more difficult and requires 
the use of semiclassical approximations [28l |. 

Our calculation of C makes use of its three crucial properties. First, C commutes with 
the space-time reflection operator VT , 

[C,VT]=0, (28) 

although C does not commute with V or T separately. Second, the square of C is the identity, 

= 1, (29) 

which allows us to interpret C as a reflection operator. Third, C commutes with H, 

[C,H] = 0, (30) 

and thus is time independent. To summarize, C is a time-independent PT-symmetric re- 
flection operator. 

The procedure for calculating C begins by introducing a general operator representation 
for C of the form U£ 



C = e^^^'P^V, (31) 

where V is the parity operator and Q{x,p) is a real function of the dynamical variables x 
and p. This representation conveniently incorporates the three requirements ()28|) - (jHUjl . 

The representation C = e'^V is general. Let us illustrate this simple representation for C 
in two elementary cases: First, consider the shifted harmonic oscillator H = jp"^ + ^x"^ + iex. 
This Hamiltonian has an unbroken VT symmetry for all real e. Its eigenvalues En = 
n + ^ + are all real. The C operator for this theory is given exactly by C = e'^V, where 
Q = —ep. Note that in the limit e 0, where the Hamiltonian becomes Hermitian, C 
becomes identical with V. 

As a second example, consider the non-Hermitian 2x2 matrix Hamiltonian ()22|1 . The 
C operator in (j24j) can be easily rewritten in the form C = e'^V, where Q = In ( |~^!|^° ) . 

Here, a2 = q*)- Again, observe that in the limit ^ ^ 0, where the Hamiltonian becomes 
Hermitian, the C operator becomes identical with V. 

We will now calculate C directly from its operator representation (jHTj) and we will show 
that Q{x,p) can be found by solving elementary operator equations. To find the operator 
equations satisfied by Q we substitute C = e'^V into the three equations - (^01) in turn. 
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First, we substitute (|HT|) into the condition (j2Hl) to obtain 

gQ(£,p) _ 'pq-QQ(x,p)'pq- _ gQ(-£,p)^ 

from which we conclude that Q{x,p) is an even function of x. Second, we substitute dHU) 
into the condition and find that 

gQ(x,p)-pgQ(z,p)p _ ^Q{x,p)^Qi-x-p) _ 

which imphes that Q{x,p) = —Q{—x, —p). Since we aheady know that Q{x,p) is an even 
function of x, we conclude that it is also an odd function of p. 

The remaining condition ()30|1 to be imposed is that the operator C commutes with H. 
Substituting C = c^^^'P^V into we get e^^^-P^fP, H] + [e^i^'P), H]V = 0. We can express 
the Hamiltonian H in (j27p in the form H = Hq + eHi, where Hq is the harmonic oscillator 
Hamiltonian Hq = ^p"^ + which commutes with the parity operator V, and Hi = ix^, 
which anticommutes with V. The above condition becomes 

2ee^^^'P^Hi = [e^^^'P\H]. (32) 

The operator Q{x,p) may be expanded as a series in odd powers of e: 

Q{x,p) = eQi{x,p) + e3Q3(x,p) + e'Q^ix^p) + ■ ■ ■ . (33) 

Substituting the expansion in (|33|) into the exponential e'^^^''^^ we get after some algebra a 
sequence of equations that can be solved systematically for the operator-valued functions 
Qn{x,p) {n = 1,3,5,...) subject to the symmetry constraints that ensure the conditions 
and (j2ni). The first three of these equations are 

[Ho,Qi] = —2Hi, 

[Ho,Q3] = -l[Qi,[Qi,m, 

[Ho,Q,] = i-,[Qu[QiAQiAQuHi]]]]-l{[Qu[Q3,Hi]] + [Qs,[Qi,Hi]]). (34) 

Let us solve these equations for the Hamiltonian in (jTTj) . for which Hq = ^p^ + and 
Hi = ix^. The procedure is to substitute the most general polynomial form for Qn using 
arbitrary coefficients and then to solve for these coefficients. For example, to solve the first 
of the equations in (jSH), [Ho,Qi] = —2ix^, we take as an ansatz for Qi the most general 
Hermitian cubic polynomial that is even in x and odd in p: 

Qi{x,p) = Mp^ + Nxpx, (35) 

where M and N are undetermined coefficients. The operator equation for Qi is satisfied if 
M = -| and = -2. 

It is straightforward, though somewhat tedious, to continue this process. In order to 
present the solutions for Qn{x,p) {n> 1), it is convenient to introduce the following notation: 
Let Sm,n represent the totally symmetrized sum over all terms containing m factors of p and n 
factors of x. For example, 5*0,0 = 1, 'S'o.s = xr", Si^i = \ {xp + px), 5*1,2 = | [x^p + xpx + px"^), 
and so on. The properties of the operators Sm,n are summarized in Ref. j22| 

In terms of the symmetrized operators Sm,n the first three functions Q2n+i are 

Qi = -IP' - 25i,2, 

Qs = W + 1*^3,2 + 85i,4 - 12p, 

Q, = - W - ¥^5,2 - ^53,4 - 64^1,6 + '-^p' + ^5i,2. (36) 
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This completes the calculation of C. Together, (|HT|) . (jH^ . and (jSHI) represent an explicit 
perturbative expansion of C in terms of the operators x and p, correct to order e^. 

To summarize, using the ansatz pip we can calculate the C operator to very high order 
in perturbation theory. We are able to perform this calculation because this ansatz obviates 
the necessity of calculating the wave functions (pn^x). 

IX. QUANTUM MECHANICS IN THE COMPLEX PLANE 

We have seen in Sec. II VI that the classical motion of particles described by PT-symmetric 
Hamiltonians is not confined to the real-x axis; the classical paths of such particles lie in the 
complex-x plane. Analogously, the new kinds of quantum theories discussed in this paper 
may also be viewed as extensions of ordinary quantum mechanics into the complex domain. 
This is so because, as we saw in Sec. IIIII the Schrodinger equation eigenvalue problem and 
the corresponding boundary conditions are posed in the complex-x plane. 

The idea of extending a Hermitian Hamiltonian into the complex plane was first discussed 
by Dyson,jwho argued heuristically that perturbation theory for quantum electrodynamics 
diverges jsoj. Dyson's argument involves rotating the electric charge e into the complex 
plane e ie. Applied to the anharmonic oscillator Hamiltonian 

H = Ip^ + Ix' + Igx' ig>0), (37) 

Dyson's argument would go as follows: Rotate the coupling g into the complex-^f plane 
to —g. Then the potential is no longer bounded below, so the resulting theory has no 
ground state. Thus, the energies En{g) are singular at g = and the perturbation series 
for En{g), which are series in powers of g, must therefore have a zero radius of convergence 
and must diverge for all 5^ 7^ 0. These perturbation series do indeed diverge, but there is 
flaw in Dyson's argument, and understanding this flaw is necessary to understand how a 
non-Hermitian PT-symmetric Hamiltonian can have a positive real spectrum. 
The flaw in Dyson's argument is simply that the eigenvalues of the Hamiltonian 

H = Ip^ + Ix' - \gx^ (^>0), (38) 

are undefined until the boundary conditions on the eigenfunctions are specified. These 
boundary conditions depend crucially on how this Hamiltonian with negative coupling is 
obtained. Dyson's way to obtain H in would be to substitute g = \g\e''^ into (fHTjl and 
to rotate from 6 = to 6 = it. Under this rotation, the energies En{g) become complex: 
The En{g) are real and positive when g > Q but complex when (? < 0. The PT-symmetric 
way to obtain is to take the limit 5: Q ^ 2 oi H = \p^ + \x^ + \gx^{ixY {g > 0). When 
(IHHj) is obtained by this limiting process, its spectrum is real, positive, and discrete. 

How can the Hamiltonian ()38|) possess two such astonishingly different spectra? As we saw 
in Sec. IIIH the answer lies in understanding the boundary conditions satisfied by the wave 
functions 4>n{x)- Under Dyson's rotation the eigenfunctions (pni^) vanish in the complex-x 
plane as — > 00 inside the wedges — 7r/3 < argx < and — 47r/3 < argx < —n. Under 
the VT limiting process, in which the exponent 6 ranges from to 2, 0„(a;) vanishes in the 
complex-x plane as |x| — 00 inside the wedges — 7r/3 < argx < and — vr < argx < —2tt/3. 
In the latter case the boundary conditions hold in wedges that are symmetric with respect 
to the imaginary axis (see Fig. |2)); these boundary conditions enforce the VT symmetry of 
H and are responsible for the reality of the energy spectrum. 
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Apart from the differences in the energy levels, there is another striking difference between 
the two theories corresponding to H in ()38p . Under Dyson's rotation the expectation value of 
the operator x remains zero. This is because Dyson's rotation preserves the parity symmetry 
of H in (|38p . However, under our limiting process, in which 5 ranges from to 2, this 
expectation value becomes nonzero because as soon as 5 begins to increase, parity symmetry 
is violated (and is replaced by VT symmetry. 

The nonvanishing of the expectation value of x has important physical consequences. 
We suggest in the next section that VT symmetry may be the ideal quantum field theoretic 
setting to describe the dynamics of the Higgs sector in the standard model of particle physics. 



X. PHYSICAL APPLICATIONS OF PT-SYMMETRIC QUANTUM THEORIES 

It is not known whether non-Hermitian, PT-symmetric Hamiltonians can be used to de- 
scribe experimentally observable phenomena. However, non-Hermitian Hamiltonians have 
already been used to describe interesting interacting systems. For example, Wu showed that 
the ground state of a Bose system of hard spheres is described by a non-Hermitian Hamil- 
tonian Wu found that the ground-state energy of this system is real and conjectured 
that all the energy levels were real. HoUowood showed that even though the Hamiltonian 
of a complex Toda lattice is non-Hermitian, the energy levels are real [3l|. Non-Hermitian 
Hamiltonians of the form H = + ix^ and cubic quantum field theories arise in studies of 
the Lee- Yang edge singularity and in various Reggeon field theory models 0. In each 
of these cases the fact that a non-Hermitian Hamiltonian had a real spectrum appeared 
mysterious at the time, but now the explanation is simple: In each case the non-Hermitian 
Hamiltonian is PT-symmetric. In each case the Hamiltonian is constructed so that the 
position operator x or the field operator is always multiplied by i. 

An experimental signal of a complex Hamiltonian might be found in the context of con- 
densed matter physics. Consider the complex crystal lattice whose potential is V{x) = 
i sin X. While the Hamiltonian H = + i sin x is not Hermitian, it is PT-symmetric and 
all of its energy bands are real. However, at the edge of the bands the wave function of a 
particle in such a lattice is always bosonic (27r-periodic), and unlike the case of ordinary 
crystal lattices, the wave function is never fermionic (47r-periodic) [s^. Direct observation 
of such a band structure would give unambiguous evidence of a PT-symmetric Hamiltonian. 

The quartic PT-symmetric quantum field theory that corresponds to i7 in (j?!) with = 4 
is described by the "wrong-sign" Hamiltonian density 

n = y{^,t) + \[vM^,t)f + i/iV'(x,t) - W{^,t). (39) 

This theory is remarkable because, in addition to the energy spectrum being real and pos- 
itive, the one-point Green's function (the vacuum expectation value of v?) is nonzero j33|. 
Also, the field theory is renormalizable, and in four dimensions it is asymptotically free and 
thus nontrivial [i^. Based on these features, we believe that the theory may provide a useful 
setting in which to describe the dynamics of the Higgs sector in the standard model. 

Other field theory models whose Hamiltonians are non-Hermitian and PT-symmetric 
have also been studied. For example, PT-symmetric electrodynamics is particularly inter- 
esting because it is asymptotically free (unlike ordinary electrodynamics) and because the 
direction of the Casimir force is the negative of that in ordinary electrodynamics [35I]. This 
theory is remarkable because it can determine its own coupling constant. Supersymmetric 
PT-symmetric quantum field theories have also been studied |36[. 
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How does a g^p^ theory compare with a gip'^ theory? A gif^ theory has an attractive 
force. Bound states arising as a consequence of this force can be found by using the Bethe- 
Salpeter equation. However, the g^p^ field theory is unacceptable because the spectrum is 
not bounded below. If we replace g by ig, the spectrum becomes real and positive, but 
now the force becomes repulsive and there are no bound states. The same is true for a 
two-scalar theory with interaction of the form igip'^Xi which is an acceptable model of scalar 
electrodynamics that has no analog of positronium. It would be truly remarkable if the 
repulsive force that arises in a "PT-symmetric quantum field theory having a three-point 
interaction could explain the acceleration in the expansion of the universe. 

We believe that the concept of VT symmetry as a generalization of the usual Dirac 
Hermiticity requirement in conventional quantum mechanics is physically reasonable and 
mathematically elegant. While the proposal of VT symmetry is unconventional, we urge 
the reader to keep in mind the words of Michael Faraday: "Nothing is too wonderful to be 
true, if it be consistent with the laws of nature." 
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